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We investigate exciton bound- state formation and crystallization effects in two-dimensional electron-hole bilay- 
ers. Performing unbiased path integral Monte Carlo simulations all quantum and Coulomb correlation effects 
are treated on first principles. We analyze diverse pair distribution functions in dependence on the layer sep- 
aration, particle density and hole-to-electron mass ratio and derive a schematic phase diagram for the neutral 
mass-asymmetric bilayer system. Our simulations reveal a great variety of possible phases namely an exciton 
gas, an exciton crystal, an electron-hole liquid and a hole crystal embedded in an electron gas. 
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1 Introduction 

Strongly correlated classical and quantum two-component systems are of high interest in many fields, for a re- 
cent overview see (TJ. Examples include electron-hole plasmas in semiconductors, e.g. [2-5], dense electron-ion 
plasmas (warm dense matter) and complex or dusty plasmas (6j|8). While the simple one-component plasma 
(or jellium) model of a Coulomb system has been studied in detail it disregards many fundamental properties, 
in particular the formation and phases of bound states. Much less is known about "real" two-component plas- 
mas. Although a number of studies have been performed in recent years, e.g. |9}|T2|, many questions remain. 
In particular, the phase diagram and its dependence on the mass ratio of the two components is known only 
qualitatively. 

Here, we analyze these questions by considering the simplest realization of a neutral two-component plasma: 
A bilayer system with spatially separated light negative and heavy positive charges, in the following referred to as 
electrons ("e") and holes ("h"). At low temperature and large mass anisotropy, M = rrih/m e ^> 1, the (classical) 
heavy component has a tendency to form a Wigner crystal whereas the light component forms a neutralizing 
background. Thus, a strong (weak) coupling situation is realized for the heavy (light) component, which can be 
quantified in terms of the one-component coupling parameter, T a = [e 2 /ea]/[h 2 /m a a 2 ] = a/as, where m a 
are the effective particle masses (a = e, h), e is the elementary charge, e the dielectric constant, and H denotes 
the reduced Planck constant. For each species, T a depends on m a (via the Bohr radius as = h 2 e/m a e 2 ) and 
the particle number density p ~ 1 /aP (a is the average inter-particle distance and D the dimensionality of the 
system). When e-h bound state (exciton) formation takes place, the residual (screened) interaction between these 
almost neutral complexes weakens, i.e. T is reduced, up to the point where a crystallization into a Wigner solid 
is prevented. Generally the formation of few-particle bound states is limited to low temperatures. Their correct 
theoretical description calls for a full quantum-mechanical treatment, in particular if temperature and density 
effects, as well as their stability against quantum fluctuations is addressed for low-dimensional systems. This 
provides the motivation of the present work: After giving a brief overview on existing results on e-h bilayers, we 
perform a first-principle numerical analysis of thermodynamically stable phases of a two-component Coulomb 
plasma in the bilayer geometry with spatially separated electrons and holes, and map out major regions of the 
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phase diagram. Thereby the focus is on the inter-layer and intra-layer particle correlations, leading to bound state 
formation and charge ordering, respectively. 

2 Electron-hole bilayers 

The main advantage of a bilayer system is the easy control of the coupling between the layers via the layer 
separation d. Actual realizations range from electron-ion quantum plasmas [1, 9] to electron-hole systems in 
semiconductors (2p3). In the latter spatially separated electrons and holes can be realized in double quantum 
well heterostructureT(3||5j. Recent advances in fabrication technology (lithography techniques and molecular 
beam epitaxy) allow for a good control of the electron/hole density and the spatial inter- well distance. Both 
classical and quantum electron-hole bilayers have been a subject of active research during the last decades [10, 
|TT][T4] - [T6] |. Classical bilayers (corresponding to high temperature and low density) have been investigated by 
molecular dynamics simulations to reveal their static and dynamical properties (15][T6) . The phase diagram has 
been explored as a function of the layer separation d = d/a, with the principle result that for 1 < d < 2.5 
the inter-layer correlations stabilize the Wigner crystals in both layers (the melting temperature is increased) 
compared to a single layer Wigner crystal. For quantum bilayers, at zero temperature, the static and dynamic 
dielectric response has been analyzed in Refs. p9f[2T| . The presence of a second (oppositely charged) layer 
was found to enlarge the region where a crystalline phase exists in the d — r s plane to higher densities, p c (d) > 
p c (d — » oo ), or r s {d) < r s (oo) (r s = Y e denotes the Brueckner parameter). The critical value of the parameter 
r s , characterizing the critical particle density, changes from r S)C « 37 for independent layers [ 22 1 to r S)C « 20 for 
a bilayer. Quantum melting of a Coulomb crystal is caused by the zero-point motion of the particles around their 
lattice sites. Similar to the classical case, these fluctuations are suppressed by attraction to the particles residing 
in the second layer, in particluar if they also form a lattice. In this way, a stack of two Wigner lattices appears to 
be more stable then a single layer crystal. 

The second component plays also a prominant role when electron-hole bound states are formed, e.g. for 
systems of spatially indirect excitons. Then the quantum bilayer is similar to a 2D system of dipoles. The 
transition from Coulomb (dominant at d ~ 1) to dipolar (d < 1) correlations shifts the critical density of the 
solid-gas phase transition to a lower value [23 ], which scales with the dipole moment (or the layer separation) as 
p c a% = l/frr^ = d 4 /290. In the dipole approximation (d <C 1), the bilayer phase diagram should be similar 
to that for bosonic dipoles in 2D and includes superfluid and normal gases (24) as well as a dipolar solid, if the 
composite particles obey Bose statistics. A more accurate treatment that accounts for the internal structure of 
bound states and allows to release the restriction d <C 1, was recently performed, see Ref. 1 25 1 . There it was 
shown that the a realistic finite layer width restores the 1 / r-Coulomb scaling at small distances which makes 
an exciton solid unstable at high densities similar to the case of a 2D Wigner crystal. Besides macroscopic 
systems, mesoscopic bilayers with full account of Fermi statistics have been studied recently [26]. Such systems 
are of particular relevance for experiments with carriers confined by an external trap potential, e.g. electron-hole 
systems in coupled quantum dots, where both the distance between the dots and the in-plane carrier density are 
tunable parameters. In contrast to bulk systems, in a mesoscopic bilayer, a change of d leads to a sequence of 
structural transitions and shell structure formation (26) similar to complex plasmas (7). 

Another important aspect, which has not been studied in detail so far, concerns the influence of the mass- 
asymmetry M between the bilayer charge carriers. A mass ratio M ^ 1 is expected to have a significant effect 
on the phase diagram, as the zero-point energies of light electrons and heavy holes differ, which in turn will 
influence inter-particle correlations (14]|27). Depending on the particle density, the light component (electrons) 
can be in a homogeneous Fermi gas phase, while the heavier component (holes) can form a Wigner lattice when 
the mass ratio exceeds a critical value, M > M c (with M c > 4 for typical electron densities in semiconductors 
which is about a factor 20 lower than in a 3D plasma |9|). The critical value M c can be tuned by the layer 
separation d (T4| . This is in favor of Wigner crystallization (at least for a heavy component) in a two-component 
Coulomb system. In contrast, in a 2D one-component system (electron gas) crystalliztation (which is expected to 
appear for r s > 37 [22]) could not yet be observed Q since the typical particle densities are too high (r s < 10), 
and quantum fluctuations destroy any spatial ordering. 



A recent experimental demonstration of Wigner islands of electrons in a mesoscopic structure covered with a helium film deals with 
the near-classical limit 1 281 . 
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In the following, we present a more systematic analysis of macroscopic bilayers, by considering a wider range 
of densities/inter-layer separations and mass ratios. We identify different structural phases by analyzing the pair 
distribution functions and reconstruct the phase diagram. For this purpose we adopt the path integral Monte Carlo 
(PIMC) technique for simulations of spatially separated electron-hole systems (14). We note that quantum effects 
between electrons and holes are fully taken into account by the PIMC formalism based on the stochastic sampling 
of the TV-particle density matrix. 

3 Theoretical model 

A wide class of electron-hole bilayers can be described by the Hamiltonian 

N e N h 2 



, 2ra a ^ e|r. 

= 1 i<3 



1,01 x 3,cx\ 



The last term in Eq. ([I]) mimics the coupling between the electrons and holes residing in different layers. Ob- 
viously, the ground-state and thermodynamic properties of the system depend on the following experimentally 
accessible and managable parameters: (i) the hole-to-electron mass ratio M = mh/m e , (ii) the inter-layer sep- 
aration d = d/dB, (iii) the in-plane particle density pa 2 B = and (iv) the temperature ksT (measured in 
Hartree: Ha = h 2 /m Q aB). In the simulations we considered a bilayer system with N Q = = 64 electrons 
and holes. The particle species are confined to their respective two-dimensional layer of zero width, separated 
by d. To simulate an infinite system we use a rectangular simulation cell with periodic boundary conditions. 
The dimensions of the cell with L y = \^3L x /2 and the number of particles were chosen to fit the symmetry 
of a hexagonal lattice, once a Wigner crystal is formed. The temperature was kept fixed to /c^T = 1/3000 Ha. 
The other parameters (M, r Sl d) were varied, each in the range 1 ... 100. More details regarding the setup of the 
PIMC simulations can be found in Ref. |T4| . In order to identify different structural phases we exploit the pair 
distribution function between different species (a/3 = ee, hh, eh) 

9ap(r) = C a(3 J2 ( S (\ r *><* - T oA ~ r ))> ( 3 ) 

* = 1 3=3af3 + l 

2 . 1 

Cee(hh) = — j— rv, j e e(hh) = h Ceh = , j e h = . (4) 

^ V e(h)UV e ( h ) - lj iV e iV h 

In particular, the formation of a Wigner solid in the electron (hole) layer is ascribed to a jump in the ratio of the 
first maximum and the next minimum in g e e(hh)> i- e - 

A = max[# aa (r)] / min [g aa (r)} , (5) 

r<2r c r<2r c 

where r c = max[^ aa (r)] is the radius of the first correlation shell. This criterion has been found sensitive to 
a solid-liquid/gas phase transition [29, 30 ]. Furthermore, to identify phases composed of excitons we use the 
fraction of bound electron-hole states, which can be estimated by the inter-layer pair distribution function g e h (for 
more details see (T7j[l8)), 



/ r[g Qh -l}< 
Jo 



X — iV eh ^ JO (( :\ 







Here r& is the second zero of r[g e h(r) — 1] on the right-hand side of the exciton peak located near r = 0. 
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4 Results and discussion 

4. 1 Characterization of different bilayer phases 

Fig. [T] presents our simulation data for the pair distribution functions (left panels) and corresponding density 
"snapshots" (right panels) for characteristic configurations of the electron-hole bilayer system. 

Classical exciton gas. Fig. [T^ exemplifies the excitonic gas phase as realized for a mass-asymmetric bilayer 
with M = 10 and an inter-layer distance much smaller compared to the intra-layer particle distance (d/r s = 
0.1 <C 1). Therefore we are in a regime where the dipole approximation is valid. Indeed, g^(r) is sharply 
peaked at r = reflecting the formation of bound states (indirect excitons). A comparison of the snapshots of the 
(instantaneous) particle density in each layer shows that all particles are bound pairwise (cf. the uppermost two 
right-hand panels). The structureless form of g e e(hh) shows that the excitons interact only weakly, i.e., they move 
almost freely and form a gaseous phase. Expressed in terms of the dipole coupling parameter D (see Ref. [ 24 ]), 
this regime is characterized by D = d 2 (M + l)/y/irr s ~ 6.2, and we expect a dipolar gas with moderate 
correlations and a superfluid transition temperature T c [Eq] = ksTc/Eo ~ 1.25 (cf. the D-dependence of T c in 
Ref. pjj scaled in units Eq = h 2 p/(M + 1)). Since T c is by factor 10 -2 smaller than the temperature simulated 
in FigTjl^, where T[Eq] = ksT/Eo — (/c^T/Ha^M + 1)tt^ 2 ~ 115, the exciton gas is essentailly classical, 
and effects of Bose/Fermi statistics are irrelevant. The situation does not change much if M = 1, because the 
degeneracy of the exciton gas remains low. 

Degenerate exciton gas/liquid. We now consider the mass -symmetric case and increase the particle density 
(reduce r s ), see Fig.JTJ). Since for d = 10 and r s = 20 the condition d/r s < 1 still holds, the excitons still behave 
as dipoles. This is supported by an analysis of the microscopic effective exciton-exciton interaction potential (25). 
For the parameters used in Fig.[lJ>, we have D ~ 5.6 and T c [Eq] ~ 1.3. The simulated temperature, T[Eq] ~ 
0.84, is now below T c . Accordingly the exciton gas is superfluid. In the quasi-2D bilayer excitonic system 
the superfluid transition will be of the Berezinskii-Kosterlitz-Thouless type pT][32| and is associated with the 
emergence of topological order, i.e. the pairing of vortices with opposite circulation. The correlations of the order 
parameter (for example, the one-particle density matrix) decay algebraically in space. A true long-range order is 
recovered only at T = 0. This scenario is applicable to a wide variety of quasi two-dimensional systems. Recent 
experimental realizations of degenerate dipolar atomic gases (33), bosonic molecules (34), and indirect excitons 
in quantum wells (3}|5) have stimulated theoretical analyses of their quantum many -body properties. The present 
bilayer system provides an interesting test case. In our system, the onset of degeneracy and spatial coherence in 
both layers (which leads to the superfluid response) is reflected in the behavior of the pair distribution functions. 
The incipient overlap of the excitonic wavefunctions (around r ~ r s ) can be noted from the particle density 
plots (right-hand panels of Fig.[lJ>). Simultaneously, we observe the build-up of short-range correlations which 
is indicated by the few (weak) oscillations in g e e(hh)- The density snapshots point to a local hexagonal ordering: 
Though noticeable spatial correlations exist between nearest neighbor and even next-nearest neighbor excitons, 
true long-range order does not appear. Such behavior is typical for correlated liquids with repulsive interactions. 
Finally, we note, that the correlations not only influence the static structural properties, but they significantly 
modify also the quasi-particle excitation spectrum (leading to the formation of a phonon-maxon-roton branch). 
The latter has a strong effect on the superfluid transition temperature (and stability of the superfluid phase to 
thermal fluctuations), in particular, in 2D dipolar systems (24). 

Excitonic crystal. An interesting situation appears when the inter-layer and intra-layer correlations are of the 
same order. If d ~ r s , the bilayer system does not behave as a pure dipolar system. For the parameters used in 
Fig. [I]: we get D ~ 124 and T[Eq] = 4.6, i.e., we are in the classical domain of the phase diagram determined in 
Ref. (24), and far away from the quantum crystallization transition observed at D ~ 17 (23) . In this case, a crude 
estimate of the melting temperature of a dipolar crystal follows from the classical dipole coupling parameter 
~ 60 [12], yielding ksT /Ha = d 2 /(T^rf) ~ 8 • 10 -4 |d =rs= 2o- This means the simulated temperature 
(fc#T = 1/3000 Ha) is 2.5 times below the melting temperature and, indeed, we observe an almost perfect 
hexagonal lattice in both layers. Obviously the bilayer geometry plays a crucial role, and we directly observe 
the stabilizing effect of the second layer mentioned in the introduction: A single electron layer would form a 
Wigner crystal only at r s > 37 (22), whereas Wigner crystal formation in a bilayer sets in at significantly higher 
densities, e.g. at about r s « 20 for the mass-symmetric case M = 1, see Refs. [19-21). In qualitative agreement, 



our data for a bilayer with heavy holes yield an exciton crystal already for r s > 10. This state is characterized 
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Fig. 1 Left column: Hole-hole (#hh, black lines), electron-electron (g ee , red lines) and electron-hole (g e h, green lines) pair 
distribution functions (averaged over 100000 Monte Carlo steps). Middle/right columns: Snapshots of the instantaneous 
particle density (averaged over 10000 MC steps) in the electron and hole layer, respectively. Varied parameters are mass ratio 
M, inter-layer distance d, and the Brueckner parameter r s . The three upper rows (from top to bottom) typify the excitonic 
phase: a) a classical exciton gas, b) a degenerate exciton gas and c) an exciton crystal. The last two rows (no e-h bound states) 
mark d) a hole Wigner crystal and e) an electron-hole plasma. All PIMC simulations were performed with N e — — 64 
particles in each layer. The temperature was kept fixed at ksT — 1/3000 Ha. 



Copyright line will be provided by the publisher 



6 



J. Schleede, A. Filinov, M. Bonitz, and H. Fehske: Phase diagram of bilayer electron-hole plasmas 



by strongly oscillating pair distribution functions in both layers. The observed difference in the peak heights of 
#hh and #ee results from the mass-asymmetry. Since the lighter electrons have a larger zero-point energy, they 
are not as localized as the holes around the hexagonal lattice sites. Furthermore we can test the appearance of 
the present bilayer exciton solid against the predictions of the exciton approximation (25) (here the analysis has 
been performed for a mass ratio M = 2.46, which corresponds to the paramaters of ZnSe-based quantum wells). 
The excitonic phase diagram derived in Ref . (25) predicts that such kind of exciton solid can exist only in a finite 
density interval, which depends on the exciton dipole moment (and on the inter-layer separation, in case of a 
bilayer). For d = 20a j g, the density simulated in Fig. [T]: (r s = 20) fits in the exciton solid density window, 
9.4 < r s < 46, estimated from Fig. 7 in (25). 

Hole crystal We now turn to the case where electrons and holes are weakly correlated. This happens at high 
in-plane particle densities (and large layer distances) when the Coulomb attraction between electrons and holes 
is screened to a large extent. Then the correlation energy within the layer, E c ~ l/r s , overcomes the inter-layer 
binding energy, Eb ~ As a result the electron-hole pairs dissociate and the spatially separated electron and 
hole subsystems undergo structural phase transitions essentially independently. For the strongly mass-asymmetric 
case, zero-point fluctuations destroy any ordering in the layer containing the light electrons, while the heavy holes 
can form a Wigner lattice, see the results for M = 100 in Fig.[T]l In the hole layer, the one-component quantum 
parameter scales with the particle density as r h s = mh/m Q r Q s , leading to rj?| r e s =i = 100 well above r s , c « 37 (22) . 
Formation of the hole crystal does not necessarily go along with a significant modulation of the electron density. 
This can be seen from the density snapshot but also from the behavior of the e-h pair distribution function, which 
stays practically constant (g e h ~ 1) for all distances. 

Electron-hole plasma. The hole crystal melts, of course, if we lower the mass of the holes. For the mass- 
symmetric case shown in Fig.[lJ, a structureless electron (hole) plasma emerges in both layers. 

4.2 Phase diagram 

After having qualitatively characterized the various phases found in the electron-hole bilayer system, we now 
construct the phase diagram in the d — r s plane, see Fig. [2] The identification of the different states and config- 
urations is based on the max-min ratio A and the bound state fraction X given by Eq. ^ and ([6]), respectively. 
The Wigner crystal phase is found above a critical value A c = 2.75 and the excitonic phase was specified for 
X > X c = 0.33. Of course, these specified critical values yield only an estimate of the exact phase boundary in 
the thermodynamic limit. The first general observation from the phase diagram is that, independent of the mass 
ratio M, excitonic bound states exist only for r s > d and r s > 5. These states are marked in Fig. [5] by filled 
symbols. A similar observation was reported in Ref. (35) . 

Large mass-asymmetry. For M = 100 and r s > d we observe some of the exciton phases discussed in 
the previous section (cf. the left most column in Fig. [2]). In the neighborhood of r s = d a bilayer exciton 
crystal is formed (filled squares). By lowering the density (r s = 100) the exciton crystal vanishes, as predicted 
by the dipole approximation valid in this regime. Once the inter-particle spacing a is increased, the dipolar 
potential energy V ~ d 2 /a 3 will decay faster than the zero-point kinetic energy K ~ h 2 /ma 2 . Hence, at some 
point quantum fluctuations will destroy any spatial ordering and initiate a phase transition to a non-degenerate 
(classical) exciton gas (filled circles), far away from the diagonal r s — d (recall that large r s implies a dilute, 
weakly interacting system). Below the diagonal (r s < d), the excitons dissociate. As a result the two layers 
are weakly coupled and exhibit essentially independent structural phases. For the scanned range of r s -values, 
1 < r s < 100, the Wigner crystal in the hole layer is the thermodynamically stable phase, due to the large 
hole mass. By contrast, the electron layer at high densities (small r s ) forms a degenerate gaseous phase up to 
r s = 5. Here an electron Wigner crystal is not observed, because the two points we have simulated (r s = 20 
and r s = 100) are either below or well above the critical density r s ^ c w 37 (22). For r s = 100, the simulated 
temperature fc^T = 1/3000 Ha is not low enough to see Coulomb crystallization of a classical electron gas. 
The corresponding classical Coulomb coupling parameter T = (e 2 /a)/&BT = (l/r s )/(fc j eT/Ha) = 30 is well 
below the critical value T c « 137, required for crystallization of a classical electron gas. 

Moderate mass -asymmetry. If we decrease the mass ratio to M — 10 the stability of the excitonic phase is not 
much affected (cf. the middle panels of Fig. [2]). For r s > d and r s > 5 their phase-space domain remains almost 
the same. At r s = d = 20 a bilayer exciton crystal forms (see also Fig.[lJ). Here we note the stabilizing effect 
of the lattice in the hole layer. Due to the coupling of the electrons to the holes arranged in a regular pattern, 
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Fig. 2 Phase diagram of an electron-hole bilayer in the d — r s plane for mass ratios M = 100 (left column), M = 10 
(middle), and M = 1 (right). The upper (lower) row gives data for the electron (hole) layer. Symbols denote different ther- 
modynamic phases: crystal (squares), degenerate gas (triangles) and classical gas (circles). Bound states (indirect excitons) 
are marked by filled symbols. The superfluid phase for 2D dipole systems is located just below the blue curve (superfluid- 
normal phase boundary) and left of the green line (superfluid-dipole crystal boundary, corresponding to D = 18) as taken from 
Ref. (24). The degenerate exciton gas, shown as filled triangles, is in the superfluid state. The melting transition boundary of 
the classical dipole phase is mareked by the red line. The dotted line is a guide for the eye showing the d — r s diagonal. 



the electron Wigner crystal is observed significantly below the corresponding critical value r SjC « 37 (above the 
density) of a single layer. Obviously, the hole lattice produces the same pinning effect as impurities studied in 
Ref. (22). In the latter case, in-plane impurities enhance the electron Wigner crystallization, i.e., they shift the 
transition density to r s ~ 7.5. The same scenario applies to bilayers with large inter-layer separation, provided 
the screening effect produced by the electrons on the hole-hole interaction can be neglected. Our simulation show 
that at higher electron densities (r s < 10) zero-point fluctuations may overcome the pinning effect of the hole 
lattice and the electron Wigner crystal melts on the diagonal d = r s (see middle column of Fig.[2j d = r s = 10). 
At an inter-layer separation d = 5 the melting of the hole lattice is accompanied by the formation of excitonic 
bound states. That is to say we observe a destabilization of the hole Wigner lattice by the electron layer due to 
screening effects. 

Mass-symmetric bilayer. In the case M = 1, obviously, all phase transitions are identical in both layers. For 
r s < 10 and r s < d we find an electron-hole plasma, otherwise a (non-degenerate or degenerate) excitonic 
gas phase, see right panels of Fig. [2] Note, that for the layer separation d = 1, exciton bound states are stable 
for r s > 10, in agreement with zero-temperature diffusion Monte Carlo results [ 10 1 . We can not confirm, 
however, the two-component plasma phase boundary reported in Ref. |T0| , simply because the corresponding 
region r s ,d < 1 was not accessible by our PIMC simulations so far. As noted before, the simulated temperature 
is not sufficiently low to observe a single-layer Wigner crystal at d ^> 1 and r s > 37. A similar argument applies 
to the exciton crystal, which should exist for layer separations above a critical value d > d c , where d c « 6a b for 
M = 2.46 (25). In the vicinity of the excitonic Mott transition more refined simulations, including an accurate 
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analysis of the fraction of ionized and bound states and a full account of the Fermi statistics (36) of electrons and 
holes would be required. 

5 Summary and concluding remarks 

In this report, we analyzed the structure and stability of possible phases of a neutral electron-hole bilayer in a 
wide range of the system's control parameters: layer separation, particle density, and mass ratio at a fixed low 
temperature. To this end we performed large-scale unbiased path integral Monte Carlo simulations. The proper- 
ties of the observed exciton gas, exciton crystal, hole crystal and electron-hole plasma phases are characterized 
by pair distribution functions and snapshots of the instantaneous particle densities. For a representative set of 
mass ratios, M = 1, 10, and 100, the phase diagram was mapped out. It shows the complex interplay between 
the intra-layer and inter-layer correlations as well as stabilization and destabilization effects of the second layer 
on the Wigner crystal in a 2D Coulomb system. Currently, extensive simulations are underway that use a denser 
grid of simulation parameters that will allow to draw the phase boundaries with higher accuracy. 

Our results are important for strongly correlated two-component charged particle systems in general (TJ and 
for complex plasmas in particular (6j|7). Common to all these systems are collective phenomena resulting in a 
rich variety of phases. Of particular interest is the possibility to control the interparticle interaction in e-h bilayers 
by changing external parameters (inter-layer separation, density). These concepts are also of interest for complex 
plasmas and for charged particles in contact with surfaces (37) where similar pinning and crystal stabilization 
effects occur due to the presence of surface defects (38). 
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